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Abstract 



Some features of a program eett6f for top quark pair production and decay into six fermions 
at linear colliders are discussed. Lowest order standard model predictions for cross sections of 
some six fermion channels and for the top quark decay width are confronted with the predic- 
tions obtained within a model with an anomalous Wtb coupling. The question of wether non 
doubly resonant background can easily be reduced by imposing kinematical cuts is addressed. 
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1 INTRODUCTION 



As a top quark is the heaviest matter particle ever observed, with mass close to the energy scale 
of the electroweak symmetry breaking, investigation of its physical properties may give hints 
towards understanding physics beyond the Stnadard Model at higher energy scales. Therefore 
precise measurements, at a level of a few per mille, of the top quark properties and interactions 
will most certainly belong to the research program of any future e + e~ linear collider pQ. In 
order to match that high precision level of the measurements, theoretical predictions should 
include radiative corrections. 

As it has been shown in j2] and will be illustrated later in this lecture, effects caused by the 
off-shellness of the it pair and off-resonance contributions may be important too, especially 
for measurements at the centre of mass system (c.m.s.) energies much above the ti threshold. 
As the t- and t-quark of the reaction 

e + e~ -> it (1) 

almost immediately decay, predominantly into bW + and bW~ , with the boson decaying 
into a fermion pair, what one actually observes are 6 fermion reactions of the form 

e + e- -> bfj&hfi, (2) 

where fi = u e , u T , u, c, fi = e~ , \T ,r~ ,d, s and /{, f 2 are the corresponding weak isospin 
partners, f[ = e~,fj,~,r~,d,s, f 2 = u e , z/ M , u T , u, c. For example, a pure electroweak (EW) 
reaction 

e + e~ -> bv^bpTv^ (3) 

in unitary gauge receives contributions from 452 Feynman diagrams, neglecting the Higgs bo- 
son couplings to fermions lighter than b, while there are only two 'signal' diagrams contributing 
to reaction (J2J) in the double resonance approximation 

e+ e - ^ t*t* -> bfj{bf 2 fi (4) 

at the same time. 

In the present lecture, some features of the updated version of a computer program eett6f 
j2] that allows for computing cross sections of reactions (J2J) to the lowest order of the standard 
model (SM), with a complete set of the Feynman diagrams, will be discussed. The question 
of wether the non doubly resonant background can easily be reduced by imposing kinematical 
cuts will be addressed. Moreover, lowest order SM predictions for some six fermion channels 
of (J2J) and for the 3 particle top quark decay width will be confronted with the predictions 
obtained within a model with an anomalous Wtb coupling. 



2 A PROGRAM 

A computer program eett6f for calculating cross sections of 6 fermion reactions (J2J relevant 
for a tt-pair production and decay at c.m.s. energies typical for linear colliders has been 
written in FORTRAN 90 The program consists of 50 files including a makefile, all stored 
in one working directory. The user should specify the physical input parameters in a module 



inprms .f and select a number of options in the main program csee6f . f . The options allow, 
among other, for calculation of the cross sections while switching on and off different subsets 
of the Feynman diagrams. It is also possible to calculate cross sections in two different narrow 
width approximations, for tt-quarks, or W /± -bosons. The program allows for taking into 
account both the EW and QCD lowest order contributions. Version 1 of the program allows 
for calculating both the total and differential cross sections of (J2J) at tree level of SM. Some 
anomalous effects beyond the SM that will be discussed later have been already implemented 
in the program. The program can be used as the Monte Carlo (MC) generator of unweighted 
events as well. 

Matrix elements are calculated with the helicity amplitude method described in 0j. Phase 
space integrations are performed with the MC method. The most relevant peaks of the 
matrix element squared, related to the Breit-Wigner shape of the W, Z, Higgs and top quark 
resonances, and to the exchange of a massless photon, or gluon have to be mapped away. As 
it is not possible to find out a single parametrization of the 14-dimensional phase space which 
would allow to cover the whole resonance structure of the integrand, the program utilizes a 
multichannel MC approach. 

Constant widths of unstable particles: the massive electroweak vector bosons, the Higgs boson 
and the top quark, are introduced through the complex mass parameters 

Ml = m 2 v - im v T v , V = W,Z, (5) 

m h = m H - im H^H, M t = m t - iT t /2, 

which replace masses in the corresponding propagators, 



4 + M t 



Afto) = — (6) 



A f(<1) = -5 7T2", S F (q) 



q 2 -M 2 H ' PKHJ q 2 ~Ml 

Propagators of a photon and a gluon are taken in the Feynman gauge. 

The EW mixing parameter may be defined either real or complex, 



■ 2 n 1 m W -In i My/ , , 

sin &w — 1 or sin W = 1 —pr. (7) 

m z 

Performing substitution © both in the s- and t-channel Feynman diagrams and taking the 
complex electroweak mixing parameter of Eq. (JJJ) leads to fulfilment of the Ward identities. As 
light fermion masses are not neglected, cross sections can be calculated without any kinematical 
cuts. 

A number of checks of the program eett6f have been performed. The reader is referred to [3] 
for details. 

In the present version of the program, an anomalous Wtb coupling has been implemented. The 
model describing departures of the Wtb interaction from the SM, caused by new fundamental 
interactions at high energies, can be specified in terms of the low energy effective lagrangian 
containing terms of dimension 6 j^j 

9 
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+ h.c. (8) 



In Eq. dHJ, P± are chirality projectors given by P± — (1 ± 75) /2, = z (7^7" — 7^7^) /2, 
the anomalous couplings and /.f are assumed to be real, and W~ should be regarded as 
an effective vector field. The SM Wtb coupling is reproduced with /+ = f% = / 2 + = an d fi 
equal to the real Cabibbo-Kobayashi-Maskawa matrix element Vtb- According to the present 
experimental data, Vtb is very close to 1 6] and /{*" is strongly constrained by the CLEO 
b — > 57 data which give /{*" ~ [7]. The couplings / 2 ~ an d /2 + are & t present only weakly 
constrained, see e.g. jH]. 



In the present section, a sample of numerical results obtained with the current version of 
eett6f will be presented. The physical input parameters that are used are the same as in |2j. 

The energy dependence of the total cross sections of Q is shown in Figure 1. The full lowest 
order cross section a is plotted with the solid line, the signal cross section ai*t* with the dotted 
line, and the cross section in the narrow top width approximation, aiu with the dashed line. 
A comparison of the solid and dotted lines shows the effects of the off resonance background 
contributions to reaction ©, and a comparison of the dotted and dashed lines shows the effect 
of the off-mass-shell production of the tt-pair. Both effects are substantial, especially at higher 
energies. Therefore they should be taken into account in the future analysis of data. 



3 NUMERICAL RESULTS 
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Figure 1: Total cross sections of e + e 



6z/ M /i + 6/i as functions of the c.m.s. energy. 



The question arises, whether those big off resonance effects and background contributions 
can easily be reduced by imposing cuts. This issue will be preliminarily addressed below for 
one specific channel of reaction ((21) by imposing the following cuts that are usually used in 



literature 



#(/,beam)> 5°, 9(q, beam) > 5°, 

8(1,1') > 5°, 0(l,q)> 5°, 

Ei > 10 GeV, E q > 10 GeV, 

m(q,q')> 10 GeV, (9) 

where q, I denote quark and charged lepton, respectively. 

Results for three differential cross sections: invariant mass distribution of a Mw-quark triple 
and energy distributions of a 6-quark and /i + of e + e" — > bv^+bdu at the c.m.s. energy of 
500 GeV with cuts Q have been compared against the corresponding results without cuts in 
Figs. 2-4. The solid histograms in Figs. 2-4 show the lowest order SM results obtained with 
the complete set of the Feynman diagrams, the dotted histograms show the contribution of 
the two ti signal diagrams and the dashed histograms have been obtained in the narrow top 
width approximation. The dotted histograms in Fig. 2 cannot be distinguished from the solid 
ones. 
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Figure 2: Invariant mass distributions of a Mw-quark triple at yfs = 500 GeV without cuts 
(left) and with cuts (right). 



By comparing the plots on the left and right hand side of Figs. 2-4, one sees that, with cuts 
(|9*]l. the background and off resonance contributions are reduced to a similar extent as the 
top pair production signal itself. Of course, this somewhat naive analysis does not preclude 
existence of a set of cuts that would keep intensity of the signal and reduced the background 
at the same time, but finding this suitable set need not be a completely simple task at all. 

eet6f gives also a possibility of looking at effects of the anomalous Wtb coupling defined by 
Eq. (jHJ that have been implemented in the program. Figs. 5 and 6 illustrate how the possible 
existence of an anomalous coupling f£ may change angular and energy distributions of t- 
quark decay products in the rest frame of the top. The solid histograms show the SM results, 
while the dotted histograms represent results in presence of the anomalous Wtb coupling with 
fi = 1, fi = j'2 = and /r; = 0.1, both in Fig. 5 and 6. In Fig. 5, the well known fact 
from the literature that angular distribution of the charged lepton resulting from the decay 
t — > bfi + u^ is the most efficient analyzer of the top-quar spin jHj is recollected. Solid histograms 
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Figure 3: Energy distributions of a 5-quark at y/s = 500 GeV without cuts and with cuts. 

in Fig. 5 clearly show proportionality of the angular distribution of [i + to (1 + cos#) if the 
spin of the decaying top-quark points in the positive direction of the quantization axis (spin 
up) and to (1 — cos#) if the spin of the decaying top-quark points in the negative direction of 
the quantization axis (spin down). 

The angular and energy distributions of a 6-quark and /i + resulting from the decay of t-quark 
produced in e + e~ annihilation at y/s = 360 GeV are plotted in Figs. 7 and 8. Again the solid 
histograms show SM results while the dotted, dashed and dashed-dotted histograms represent 
results in presence of the anomalous Wth coupling. One sees that the flat histograms of 
Fig. 5 corresponding to decay of the unpolarized top quark have developed a slope due to 
the Lorentz boost from the rest frame of t-quark to the c.m.s. of e + e~. Similarly, how the 
Lorentz boost changes the energy distributions of the top decay products is illustrated in 
Fig. 8. It is interesting to note, how the relatively big effects arising if only a single anomalous 
coupling is nonzero, i.e. f£ = 0.1 and / 2 ~ — (dotted histograms), or / 2 + = and f% = 0.1 
(dashed histograms), are reduced by interference, if f£ = 0.1 and / 2 ~ = 0.1 (dashed-dotted 
histograms) . 



4 SUMMARY AND OUTLOOK 



Top quark pair production and decay into 6 fermions in e + e~ annihilation at c.m.s. energies 
typical for linear colliders can be theoretically studied to the lowest order of SM with a program 
eett6f . A sample of results on selected channels of reaction (J2J) have been presented. It has 
been shown that, although the two ti signal diagrams dominate total cross sections even at 
c.m.s. energies much above the ti threshold, the effects related to off-mass- shell production 
of the tt-pair and the off resonance background may be relevant for the analysis of future 
precision data. It has been illustrated that these effects are not reduced by imposing typical 
cuts given by 0. Some extensions of the SM have been implemented in the program and 
results illustrating effects of the nonstandard Wth coupling Eq. (JHJ) have been shown. 
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Figure 4: Energy distributions of a /i + at = 500 GeV without cuts and with cuts. 



Work towards completing implementation of all the channels of (J2J) possible in the SM is 
underway. In order to match high precision of future experiments it is also required to im- 
plement higher order effects. As it is not feasible to calculate radiative corrections to the full 
set of Feynman diagrams that contribute to any specific channel of (J2J), it would be desirable 
to include higher order effects at least to the two signal diagrams. This approach would be 
justified by the fact that the total cross section of reactions (J2J) is dominated by the doubly 
resonant signal. Work in this direction is planned in collaboration with the Zeuthen-Bielefeld 
group [TUj . 
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Figure 5: Angular distributions of a top quark decay products in the top rest frame. 
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Figure 6: Energy distributions of a top quark decay products in the top rest frame. 




Figure 7: Angular distributions of a 5-quark (left) and ji + (right) at y/s = 360 GeV. 
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Figure 8: Energy distributions of a 5-quark (left) and fi + (right) at y/s = 360 GeV. 
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